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The double-slit experiment strikingly demonstrates the 
wave-particle duality of quantum objects. In this famous 
experiment, particles pass one-by-one through a pair of 
slits and are detected on a distant screen. A distinct 
wave-like pattern emerges after many discrete particle im- 
pacts as if each particle is passing through both slits and 
interfering with itself. While the direct event-by-event 
buildup of this interference pattern has been observed 
for massive particles such as electrons 1 1 -3|, neutrons |4|, 
atoms |5| and molecules |6, 7|, it has not yet been mea- 
sured for massless particles like photons. Here we present 
a temporally- and spatially-resolved measurement of the 
double-slit interference pattern using single photons. We 
send single photons through a birefringent double-slit ap- 
paratus and use a linear array of single-photon detectors 
to observe the developing interference pattern. The anal- 
ysis of the buildup allows us to compare quantum me- 
chanics and the corpuscular model described in Ref. |8|, 
which aims to explain the mystery of single-particle inter- 
ference. Finally, we send one photon from an entangled 
pair through our double-slit setup and show the depen- 
dence of the resulting interference pattern on the twin pho- 
ton's measured state. Our results provide new insight into 
the dynamics of the buildup process in the double-slit ex- 
periment, and can be used as a valuable resource in quan- 
tum information applications. 

Interference pattern measurements using individual pho- 
tons have been carried out with relatively slow exposing 
charge-coupled device (CCD) cameras ll9UTT1l and by scan- 
ning a single-photon detector through a detection plane 1 12|. 
However, CCD cameras do not resolve the impact of individ- 
ual photons, and scanning single-photon detectors cannot si- 
multaneously record full spatial and temporal information. In 
our setup, we use an array of 32 single-photon avalanche de- 
tectors (SR\D) |[T3l[T3| as a detection "screen" for our double- 
slit setup. Using this SPAD array in our interference setup, we 
are able to observe the buildup of the double- slit interference 
pattern with high resolution in both space and time. 

Our experimental setup, shown in Fig. [T] uses a Sagnac- 
type source 1 15 1 producing photon pairs, a vertically -polarised 
photon at 842 nm and a horizontally-polarised photon at 
776 nm, via the nonlinear process of spontaneous paramet- 
ric downconversion (SPDC) I1I61I . The 776 nm photon passes 
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FIG. 1. Experimental setup. The Sagnac-type source produces po- 
larization entangled photon pairs. One photon is coupled into single- 
mode fiber. A birefringent calcite crystal displaces photons with hor- 
izontal polarization, and a crystal (CC) compensates for path length 
difference. A polarizer (P) erases any distinguishing information 
about the photons. Two lenses (LI and L2) determine beam size, and 
a third lens (L3) focuses the beam vertically onto the SPAD detectors. 
The other photon is sent through a polarization analyser consisting of 
a half wave plate (HWP), quarter- wave plate (QWP) and polarizing 
beamsplitter (PBS). It is then coupled into one of two single-mode 
fibres connected to detectors (DI and D2). The inset shows a photo 
of the 32-pixel SPAD array |13 14 1. 



through a series of waveplates and a polarizing beamsplitter. 
It is detected by D2 and used to herald the presence of the 
842 nm photon 1 17|. The 842 nm photon is coupled into a 
single-mode fibre, and a polarization controller (not shown) 
prepares the state in an equal superposition of horizontal (H) 
and vertical (V) polarizations. This is then outcoupled, re- 
sulting in a free- space Gaussian spatial mode with a waist of 
1.3 mm. This beam is collimated and sent to a polarization- 
based double slit composed of a calcite beam displacer. The 
birefringence of this crystal results in the displacement of hor- 
izontally polarized photons by 3.68 mm with respect to the 
vertically polarized photons. The beam displacer maps the 
polarization state of a photon into a spatial state, which is en- 
coded in its path. These two paths are analogous to a double- 
slit apparatus. They are orthogonally polarized and thus carry 
distinguishing information, which is erased by a polarizer set 
at 45 degrees. A compensating crystal (CC) is placed after 
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FIG. 2. Interference pattern buildup. Panel (A) shows first 200 her- 
alded counts in time, and panels (B-D) depict the statistics of the first 
20, 200 and 2000 heralded detections. 



the beam displacer to make the two path lengths equal, and a 
series of lenses maps the interference pattern onto the SPAD 
array. 

Each of the 32 detectors in the SPAD array records the ar- 
rival time of single photons with a timing uncertainty of about 
150 ps, which is the combined timing jitter of the detectors 
and time tagging logic. Fig. [2jA) shows the arrival times of 
the first 200 detection events passing through the slits. The ac- 
cumulation of these events results in an interference pattern, 
as shown in Fig.|2jB-D). After the detection of 2000 photons, 
the interference pattern becomes very clear, with a visibility 
of 93 ± 2%. This visibility is not perfect as a result of inexact 
compensation of the two path lengths. A movie and additional 
measurements using a coherent source can be found in Sup- 
plementary Information. 

Our ability to accurately measure the arrival times of pho- 
tons allows us to test the predictions of an alternative corpus- 
cular theory, designed to explain the phenomenon of interfer- 
ence without wave-particle duality |8|. In this theory, detec- 
tors are modelled as deterministic learning machines, which 
are able to reproduce the interference pattern after many pho- 
ton detections. The detectors' internal states update after each 
photon detection (see discussion in Supplementary Informa- 
tion), improving their knowledge of the pattern. 

Using two statistical methods and the measured buildup of 
the interference pattern, we examine the predictions of this 
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FIG. 3. Statistical tests. (A) Coefficient of determination. For a 
given photon number, the statistics of R^ is generated after 10^ nu- 
merical Monte Carlo simulations for the corpuscular and quantum 
mechanical models. The red (blue) belt shows 50% of the most fre- 
quent values of R^ for the case of the quantum mechanical (corpus- 
cular) model. The red points are the values of R^ for experimental 
data. (B) Likelihood ratio test. The smallest likelihood ratio value 
is log A = 0.83, which shows that quantum mechanics is a better 
indicator of the behaviour seen in nature. 



corpuscular theory and quantum mechanics. The coefficient 
of determination 1 18 1, i?^, allows us to evaluate how well each 
model predicts the final interference pattern with increasing 
detection number, while the likelihood ratio test (191, ^^ al- 
lows us to compare the two models. 

We begin by calculating the coefficient of determination, 
R'^, to see how quickly the measured data, the corpuscular 
model and quantum mechanics each reproduce the final inter- 
ference pattern. This pattern is derived from classical wave 
mechanics, and intensity is used as the only fit parameter 
(see analysis of interference pattern in Supplementary Infor- 
mation). Our experimental data gives us R'^ = 0.96 after 
190 ± 5 detections, as shown in Fig. [3jA). This tells us that 
the interference pattern is clearly visible after only 190 detec- 
tion, which we then use as a reference for comparison with 
the two models. Next, we use the Monte Carlo method to run 
10^ numerical simulations of 1 . . . 2000 photon detections for 
quantum mechanics and the corpuscular model. The statistics 
of R'^ for these simulations are shown in Fig.[3[A). Although 
both quantum mechanics and the corpuscular model eventu- 
ally predict the final pattern very well, they require 200 ±5 and 
1000 ± 10 photons, respectively, to achieve R'^ = 0.96. While 
it is clear that these statistics for the quantum mechanical sim- 
ulations and experimental data have similar trends, the coef- 
ficient of determination cannot conclusively say which model 
is better. 

In order to compare the two methods, we perform a like- 
lihood ratio test. This test tells us which model is better at 
reproducing the observed data (see discussion in Supplemen- 
tary Information). First, we calculate the probability distribu- 
tion of photon detections based on quantum mechanics. We 
then numerically simulate the corpuscular model 2.6 x 10^ 
times using the best algorithm in Ref. |8 1 to obtain its detec- 
tion probability distribution, which is dependent on the num- 
ber of detected photons. In contrast, the quantum mechani- 
cal distribution has no such dependence. Next, we calculate 
how likely it is that our experimental data emerges from these 
probability distributions and compare them using the likeli- 
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FIG. 4. Interference. (A) The round (square) points show the interference pattern of the first 2000 photons heralded by D2 (Dl). The 
triangular points show the envelope that results from heralding by either polarization. Limitations of electronics resulted in fewer coincidences 
at detectors 7 and 10. (B,D) Interference pattern fringes move as the phase is changed remotely by the QWP. The measurements are taken 
every 10-degree rotation. See Supplementary Table I for the visibilities of each set of measurements. (C,E) The trajectory of the Bloch vector 
related to the remotely prepared states heralded by (C) Dl and (E) D2. 



hood ratio, A. This test is independent of the numerical value 
of the ratio; as long as log A > 0, we can say that quantum 
mechanics is more accurate than the corpuscular model. Since 
log A > 0.83 for all points in Fig. |3jB), we conclude that 
quantum mechanics is a better indicator of the behaviour seen 
in nature. 

In a second experiment, we use the same setup and 
generate polarization-entangled photons in the state |V^) = 
7i(l^^)s,i + \^^)s,i) with fidelity 0.94. Here s,z repre- 
sent the signal and idler photons. The orthogonal polariza- 
tion states of the 842 nm signal photon, \H) and \V), are 
transformed into the spatial states |t) and ||) by the cal- 
cite crystal. These refer to the two possible paths through 
the beam displacer. The resulting entangled state is |V^) = 
T^di H)s,i + It y)s,i)' The 776 nm idler photon is sent to 
a polarization analyzer, which consists of waveplates, a polar- 
izing beamsplitter and two detectors (see Fig.[T]). The orienta- 
tion of the HWP is set such that detection by Dl and D2 corre- 
spond to projection on {\H) + |F))/V2 and {\H) - \V))/V2, 
respectively. 

After taking data for 60s, we filter the detection events 
by choosing detections at either Dl or D2 as the trigger. 
If we choose Dl as the trigger, we herald the state (|t) + 
\D)/^/2, which leads to the interference fringes shown in 
Fig. |4jA). Similarily, triggering by detection at D2 heralds 
(It) — II) )/a/2, resulting in a complementary interference pat- 
tern. The fringes are complementary because of the phase 
difference between the states heralded by Dl and D2. If we 
instead choose to herald using Dl 6>r D2 without distinguish- 
ing between the two, there is no interference pattern. This 
is because we effectively ignore the polarization state of the 
trigger photon, leaving the signal photon in a mixed state. 



Because the photons are entangled, the phase of the inter- 
ference pattern is correlated with the polarization state of the 
signal photon. To show that we indeed have entanglement be- 
tween spatial and polarization degrees of freedom, we rotate 
the QWP in the polarization analyser. The resulting effect on 
the fringes are shown in Fig.|4](B,C). The phase of the pat- 
tern is clearly dependent on the polarization state of the trig- 
ger photon. In contrast, the polarization state of the trigger 
photon would have no effect on the phase of the interference 
pattern if these were non-entangled pairs. This heralding can 
also work in reverse. By post-selecting on a particular point 
in the interference pattern, it is possible to prepare the idler 
photon in a specific polarization state. Such a flexible remote 
state preparation could be very helpful in photonic quantum 
information processing. 



The double-slit experiment, which is at the "heart of quan- 
tum mechanics", has played a central role in our understand- 
ing and interpretation of quantum theory 1 20] . Now, over two 
hundred years after the first experiments by Thomas Young 
1211 I22II . our results provide the most complete picture of 
single-photon interference to date. Additionally, our time- 
resolved measurement techniques will dramatically decrease 
the difficulty of directly measuring the wave function of a sys- 
tem by performing weak measurements 12311241 . It will also 
allow us to herald a variety of polarization states in a multi- 
plexed fashion, as well as facilitate the encoding and transfer 
of information using the hyper-entanglement of the spatial, 
temporal and polarization degrees of freedom 1.1 L, 25, ■26.1 . 



Acknowledgements 

The authors acknowledge funding from NSERC, Ontario 
Ministry of Research and Innovation (ERA program), CIFAR, 
Industry Canada and the CFI. The research leading to these re- 
sults has received funding from the European Union Seventh 
Framework Programme (FP7/2007-2013) under grant agree- 
ment no. 257646. The authors thank Simone Bellisai and 
Franco Zappa for their contribution in developing the detec- 
tors array and Christopher Granade for insightful discussions. 
PK acknowledges support by Mobility Plus project financed 
by Polish Ministry of Science and Higher Education. 



Methods 

Experimental setup. The details of the Sagnac-type 
source of photon pairs are described in Ref. |15|, with a few 
modifications. The pump is a 404 nm laser diode (Toptica 
Bluemode), and the down conversion crystal is a 30 mm 
PPKTP crystal phasematched to produce photons at 776 
and 842 nm. The output of CI has beam waist 1.3 mm. 
The calcite crystal is 41 mm long, and the compensation 
crystal is 5 mm long. Lens LI is plano-convex (f=150 
mm), lens L2 is aspherical (f=ll mm) and lens L3 is a 
plano-convex cylindrical (f=25 mm). Dl and D2 are Perkin 
Elmer SPCM-AQ4C single photon detectors. The photon 
source produced around 2 x 10^ photon pairs/second which 
resulted in around 36 x 10^ fiber coupled pairs/second. Then 
the transmission of the calcite system decreased this number 
to approximately 72 x 10^, which results in around 2000 
detected coincidences/second. The SPAD array detector 
dark count rate gives rise to approximately 5 accidental 
coincidences/second. All 32 channels of the SPAD array 
are recorded individually as time tags by two logic units 
(UQDevices). 

SPAD array. The SPAD array is a 32x1 array of single- 
photon avalanche diodes (131 [141, with pixel pitch of 100 /im 
and photon detection efficiency 5 % in the range 770 - 840 
nm. It has active area diameter of 50 fim and a dark count 
rate of 100 counts/s per pixel. For technical reasons, we use 
28 of the pixels. 
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Supplementary information 

Movie can be find here: 

Ihttp : //youtu . be/HllhJWIcUYO| 

Additional measurements - Coherent source 

Another experimental setup depicted in Fig. [5] has three 
components: an attenuated laser as a single photon source, 
a slit system |26| and a SPAD array detector [|13[il4il . 




FIG. 5. Attenuated coherent states. The Ti: Sapphire laser pumps the 
BBO. Photons are coupled into single-mode fiber using lens LI and 
collimated by L2 before passing through the double slits. Lens L3 
focuses the beam on the SPAD array. 

The attenuated coherent states experiment is based on fiber- 
coupled frequency-doubled 1 ps pulses at 396 nm attenuated 
such that the SPAD array detects approximately 200 pho- 
tons/s. This results in 0.06 average photons per pulse which, 
assuming Poissonian statistics, makes the probability of more 
than one photon arriving at the slit at the same time negligible. 

The Ti: Sapphire laser outputs 792 nm pulses, and is used to 
pump a 2 mm BBO crystal. Single mode fiber and collimat- 
ing lens L2 output a gaussian beam with a radius of 0.85 mm 
(FWHM). Lens L3 has a focal length of 10 cm. The slits are 
500/im high, 30/im wide and separated by 100 /im. 

Lens L2 is chosen such that the impinging photons' spatial 
mode size is much larger than the slits' characteristic size. 
This results in a uniform illumination and transmission of 4%. 

After passing through the slits, the photon propagates 
through lens L3 and is detected by a SPAD array in the fo- 
cal plane, where the interference pattern is formed. The pho- 
tons' spatial and temporal modes, in addition to the optics, are 
chosen such that the interference pattern's characteristic size 
is comparable to the dimensions of half of the array. This, in 
conjunction with the high quantum efficiency of the SPAD 
array at this wavelength, which is 40%, allows an optimal 
signal-to-noise ratio at the single-photon level. It also allows 
the expected interference pattern minima (maxima) to coin- 
cide with odd (even) numbered SPAD array pixels. 

The timing information of a detected photon is recorded 
only if it is coincidental with a reference pulse from the 
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FIG. 6. Interference pattern buildup for 2 slits, (a) First 200 detec- 
tions in time and (b-d) statistics of first 20, 200 and 2000 detections 
are presented. 



Ti: Sapphire pulsed laser, thus reducing background noise 
from dark counts and stray light. The timing resolution of 
the electronics and SPAD array result in a detection window 
of 312 ps and an effective 6 total dark counts/s. 

We observe the buildup of the interference pattern in time, 
as seen in Fig. [6] Note that odd-numbered SPAD pixels detect 
significantly fewer photons than even-numbered pixels, which 
correspond to the minima and maxima of the interference pat- 
tern. The contrast between the neighbouring pixels is clearly 
visible after only 20 photon detections, indicating the exis- 
tence of some sort of interference pattern from the beginning 
of the measurement. 

Additional measurements are done for three slits using at- 
tenuated coherent states. The imaging optics are adjusted to 
fit the SPAD array and interference pattern characteristic di- 
mensions. The consecutive detection events as a function of 
detector number and time, as well as the histograms of the 
recorded time tags, are depicted in Fig. [7] 



Interference pattern analysis 

After passing through the system of a birefringent crystal 
and polarizer, a wave is in superposition of two spatial gaus- 
sian modes displaced by the distance d. Assuming a charac- 
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FIG. 7. Interference pattern buildup for 3 slits, (a) First 200 detec- 
tions in time and (b-d) statistics of first 20, 200 and 2000 detections 
are presented. Note that due to technical problems, pixels 16 and 17 
were disconnected. 



teristic radius of w, the corresponding mode functions are: 

f-(x±d/2)^ 
u±{x^ z = 0) (X exp ' 



(1) 



Note that these modes are nearly orthogonal if the displace- 
ment d is sufficiently large. In our case, w = lA mm and 

d = 3.68 mm and the overlap is \/:^e~^ = 0.0096. 

Next, the standard Fresnel propagation allows to compute 
the filed in the focal plane of the lens: 

7Tx{7Tw'^X±2idfX^ 

u±{x,z = f) (xe f^ . (2) 

By modification of the relative amplitudes and the phases of 
the two modes, one can prepare any superposition of the fol- 
lowing form: 



a+'U+(x, z) + a-U-{x^ z). 



(3) 



Equations (2) and (3) allow us to predict the interference pat- 
ten shape and its dependence on the input state. 



The corpuscular model 

The corpuscular model of the double-slit experiment fH 
gives an alternative description for the buildup of the interfer- 
ence pattern. The detectors are based on deterministic learn- 
ing machines, whose internal states are updated with each 



"messenger" (photon) detection. The messengers propagate 
for a specific time after passing through the slits, acquiring a 
phase (j), which then updates the detectors' states according to 
the equations IISl : 

Pk 



7(1 -Wk-l), 


(4) 


l^k-iPk-1 + (1 - /i/c-i)e/c, 


(5) 


, /I ^\\Pk -Pk-l\\ 


(6) 



z 

where k, and 7 are constants associated with the detectors, pk 
is a parameter that is updated with each detection, and w^ is 
the internal state of the detector. 

In addition to the likelihood ration test presented in Fig.|3] 
we comment on few observations on the corpuscular model 
based on the numerical simulations and our measurements. 

We compare the this model to the other aspects of our ex- 
periment, including complementary fringes and the shifting 
interference patterns. It is very clearly shown in Fig.Qa) that 
the two sets of fringes are extracted from the same measure- 
ment data. This situation would confuse the detectors' learn- 
ing process, thus telling us that entanglement resides outside 
of the scope of the corpuscular model. 



Interference pattern visibilities 



QWP [deg] Vis., [%], Dl Vis.,[ %], D2 






93±2 


96±2 


10 


94±2 


98±2 


20 


93±2 


94=b2 


30 


97il=2 


92±2 


40 


94±2 


94±2 


50 


96±2 


93±2 


60 


90±2 


87±2 


70 


84il=2 


85±2 


80 


75±2 


89±2 


90 


91±2 


94±2 


100 


93±2 


95±2 



TABLE I. Visibilities for patterns obtained by QWP rotations, see 

Fig.g 



Likelihood ratio test 

The likelihood ratio test fJM allows to qualitatively evalu- 
ate which of the two models is a better predictor of the mea- 
sured data. This can be done by looking at the ratio, A of the 
probability, P ( I) | Ml), to get a certain set of data, D, under 
one model. Ml, to the probability, P{D\M2) to get the same 
data, D, under a second model, M2. It is convenient to take 
the logarithm of this ratio: 

log A = \ogP{D\Ml) - \ogP{D\M2), (7) 

where D = {/ci, /c2, ..., /c28} in our case represents the photon 
counts distribution measured by the SPAD array (ki, stands 



for counts of ith detector) and M = {pi,P2, ••-,^28} is the 
probability distribution given by a model ( Pi is the probability 
of detecting a photon by ith detector). If log A > 0, we can 
say that Ml describes experimental data better than M2 does. 
The probability to get a distribution, D, assuming probabil- 
ity distribution, M, is calculated using a multinomial expan- 
sion: 



P(D|M) = Yl 



n 



/ci+fc2 + ...+/e28 l<n<28 



Pn 



(8) 



The probability distribution for the quantum mechanics. 
Ml, is computed in the following way. We fit the photon 
statistics acquired after detection of 98000 photons to the 
quantum mechanical model. The fitted parameters are inten- 



sity, transverse shift and magnification. The remaining setup 
parameters were fixed to the measured values. The coefficient 
i?^ is 0.99. Note that within this model the probability distri- 
bution does not depend on detected photon number. 

This is the feature that differs quantum mechanics and the 
corpuscular model. To calculate the photon detection prob- 
ability distribution, M2, for the corpuscular theory, 2.6 x 
10^ simulations were made for photon number in the range 
1 . . . 2000. In contrast to quantum mechanics, the probability 
distribution for this model depends on the number of photons 
that have been detected. 

It is important to note that it is meaningless to compare the 
values of the likelihood ratio for different sets of data, e.g. 
different numbers of photons. 



